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ABSTRACT 

A very high resolution (R — A/AA = 200,000), high signal-to-noise ratio (S/N ~ 340) blue-green 
spectrum of the very metal-poor ([Fe/H] ~ —2.6) red giant star HD 122563 has been obtained by us at 
McDonald Observatory. We measure the asymmetries and core wavelengths of a set of unblended Fe I 
lines covering a wide range of line strength. Line bisectors exhibit the characteristic C-shape signature 
of surface convection (granulation) and they span from about 100 ms -1 in the strongest Fei features 
to 800 ms -1 in the weakest ones. Core wavelength shifts range from about —100 to — 900 ms -1 , 
depending on line strength. In general, larger blueshifts are observed in weaker lines, but there is 
increasing scatter with increasing residual flux. Assuming local thermodynamic equilibrium (LTE), 
we synthesize the same set of spectral lines using a state-of-the-art three-dimensional hydrodynamic 
simulation for a stellar atmosphere of fundamental parameters similar to those of HD 122563. We find 
good agreement between model predictions and observations. This allows us to infer an absolute zero- 
point for the line shifts and radial velocity. Moreover, it indicates that the structure and dynamics 
of the simulation are realistic, thus providing support to previous claims of large 3D-LTE corrections, 
based on the hydrodynamic model used here, to elemental abundances and fundamental parameters 
of very metal-poor red giant stars obtained with standard 1D-LTE spectroscopic analyses. 
Subject headings: stars: atmospheres — stars: population II — stars: individual (HD 122563) 



1. INTRODUCTION 

The red giant star HD 122563 is the brightest very 
metal-poor ([Fc/H] < —2.5) object known (V = 
6.2 mag). Its extreme meta l defic iency was first iden- 
tified by Wallerst ein et al.l ()1963D and since then a 
large number of studies concerning its detailed chem- 
ical composition have been published (e.g.. iSneden 



1973HLambert et al.lll974tlRvan et al.lll996HCavrel et al 



20041 ). Twenty entries more recent than 1990 are listed 
for H P 122563 in the PASTE L database of stellar param- 
eters (jSoubiran et al.ll2010D . which give: T cS = 4602 ± 
72K, logg = 1.25 ±0.22, and [Fe/H] = -2.64 ± 0.12 
(mean and standard deviation). 

Most elemental abundance studies use one-dimensional 
model atmospheres in hydrostatic equilibrium. The at- 
mospheres of cool stars, however, are far from static. Ve- 
locity and intensity fluctuations due to sur face convec- 
tion a re directly observed in the Sun (e.g., iNesis et all 
119921 : iRutten et afl 120041 ) and have been detected in- 
directly using hi gh resolution spectra of a number of 
other stars (e . g. . lAllende Prieto et al.1 120021 : iGravl 120051 : 
iDravinsI [20081 iRarmrez et al.l l2008t T In the solar pho- 
tosphere, bright granules containing hot ascending gas 
are observed separated by intergranular lanes of cool de- 
scending gas. The correlation between temperature and 
velocity fields, as well as the granule/intergranule con- 
trast in velocity and intensity, shape the spectral fea- 
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tures so that the line cores of disk-integrated profiles are 
shifted with respect to their re st wavelengths and t he line 
profiles are asymmetric (e.g.. iDravins et aT1 [l98lK The 
latter are granulation signatures that can be observed 
both in the solar spectrum and in the spectra of distant 
stars. 

Hydrodynamic simulations provide the most real- 
istic approach to the model i ng of surface inhomo- 
geneities (e.g. lNordlundlll982t iSteffen fe Frevtag] 119911: 
iLudwig et aDl2009t iNordlund et al.ll2009t ). although non- 
hydro dynamic modeling is also possible (e.g., IDravinsI 
1990). Three-dimensional model atmospheres for the 
Sun, Procyon, and K-dwarfs hav e been tested against 
high quality spectroscopi c data (lAsplund et al.l 120001 : 
lAllende Prieto et al.l 120021: IRarmrez et al.l I2009D . show- 
ing in general good agreement. These tests are crucial 
because they validate the models and allow us to use 
them with confidence to estimate the impact of "3D ef- 
fects" on stellar abu ndance and parameter determina- 
tions (jAsplundl 120051) . 

The most extreme abundance corrections proposed in 
connection with granulation effects occur at low metal- 
licities. This is mainly because the difference between 
the ID and 3D predictions for the mean temperature 
of the outer layers of metal- poor stars is very large , 
on the order of 1000 K (e.g.. ICollet eTaTl 12009. 120071 
I Gonzalez Hernandez et all |2010|) . This large tempera- 
ture difference is very important for the formation of 
features tha t are t emperature sensitive. For example, 
ICollet et al.l ()2009f ) report 3D corrections of ~ — ldcx 
(i.e., one order of magnitude) for the abundance of nitro- 
gen and oxygen as determined from molecular features, 
similar to the corrections required for the iron abundance 
determined from low excitation potential Fe I lines. They 
also estimate a correction of about +0.4 dex for the iron 
abundance determined from Fe n lines. 
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We can only rely on the corrections suggested by 3D 
models if the hydrodynamic simulations are realistic. 
Tests of the hydrodynamic simulations include measure- 
ments of the profiles of stellar absorption lines both with 
respect to shape and absolute wavelength. In this Let- 
ter, we present new spectroscopic observations for the 
very metal-poor red giant star HD 122563 and compare 
them to state-of-the-art 3D model predictions. 

2. SPECTROSCOPIC OBSERVATIONS, LINE DATA, AND 
MEASUREMENT OF LINE ASYMMETRIES 

We observed HD 122563 on June 21-25, 2010 with the 
Robert G. Tull coude spectrograph on the 2.7 m Har- 
lan J. Smith Telescope at McDonald Observatory. Us- 
ing a narrow s lit (0.34 arcsec), the spectral resolution is 
R = 200,000 (|Tull et al.lll995h . The wavelength cover- 
age of a single exposure consists of 16 non-overlapping 
pieces of ~ 20 A each. To increase wavelength coverage 
we rotated the spectrograph grating to three different po- 
sitions (setups) . In this way we were able to cover about 
70 % of the 4425 to 5575 A spectral window. To increase 
the signal-to-noise ratio (S/N), multiple 20- minute expo- 
sures of each setup (from 8 to 16, depending on the setup) 
were acquired. The total exposure time for all setups was 
about 10.5 hours. The S/N of the coadded spectrum 
increases from about 180 at 4425 A to 440 at 5575 A, 
with a median S/N = 340. The spectrosco p ic dat a were 
reduced and coadded as in iRarmrez et all (|200g ). The 
RMS of the wavelength calibration for each exposure is 
about 5xl0" 4 A (~ 30ms -1 ). Continuum normalization 
was done by fitting polynomials to the upper envelope of 
the data, first in the dire ction of dispers i on an d then 
perpendicular to it, as in IBarklem et all ((2002T) . The 
barycentric correction was determined with a precision 
better than 1ms" 1 using a program by S. Els (private 
communication) . 

A radial velocity correction was applied to each 
setup by comparing the observed core wavelength shift 
versus line strength relation with the theoretical one 
(c f. Sect - Ill). T his procedure is similar to that described 
m Gray (2009) but using a theoretical prediction instead 
of a solar reference. Basically, the zero point of the ob- 
served core wavelength shifts is modified until the av- 
erage difference with the model, for all lines in a given 
setup, is exactly zero. The absolute radial velocity we 
derive by applying this method to the three setups is 
V r = —25.39 ± 0.09 kins -1 (the error bar corresponds 
to the 1 a standard deviation of the three velocities de- 
rived) Q The uncertainty in this velocity is most likely 
due to instrumental and temperature/pressure instabili- 
ties (the R. G. Tull spectrograph is not in a vacuum). 

We identified 58 unblended Fei lines in our HD 122563 
spectrum. For these features we adopted the labora- 
tory wavelengths measured by iNave et ail (|1994l ) , which 
have a precision of about 20 ms" 1 , as estimated from 

1 This absolute radial velocity is not corrected for the star's 
gravitation al redshift, which is, however, expected to be relatively 
small (e.g., Dravins et al. 1999). Metal-poor red giants have typi- 
cal masses of ~ 0.75 Mq while we estimate a radius o f ~ 23 Rp, for 
HD 122563 based on its observed Hipparcos parallax ( v an Leeuwenl 
120071) a nd the angular diameter ob tained with the infrared flux 
method (Ramirez & Mclcndcz 2005). These values lead to a gravi- 
tational redshift of 633 m s" 1 x[(M/Mq)/(R/Rq)] ~ 0.021 kms" 1 
which is smaller than our uncertainty in V r . 
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Figure 1. Bisectors of 58 Fei lines in our spectrum of HD 122563. 
Filled circles show the location of the line cores. The dashed line 
is at zero velocity. 

the uncertainty in the wavenumber of the relevant en- 
ergy levels. The line synthesis presented in Sect. [3] was 
done using transition probabilities measured in the lab- 
oratory by a number of groups and homogeni z ed fo l- 
lowing the recommendations by lLambert et ail |l996). 
The van der Waals da mping constants used are those by 
IBarklem et al.l (j2000l ). except for five lines in which we 
used the values obtained from the classical Unsold for- 
mula. 

Line bisectors were measured in our HD 122563 spec- 
trum (and also in an identical fashion in the theoreti- 
cal line profiles, Sect. |3|) by first resampling the spectral 
lines, using cubic spline interpolations, to iso-flux points 
on both wings of the lines. The midpoints of horizontal 
segments across the spectral line wings, which define the 
line bisector, were then computed. The wavelengths of 
observed spectral line cores were determined using the 
7 data points closest to the flux minimum, which corre- 
sponds to a window of about 3.8 kms" 1 , and fitting a 
fourth order polynomial. Core wavelength shifts (here- 
after given as a velocity, i.e., Av c = cAA/A, where c 
is the speed of light and A A = A Q b s — Ai a t>) were then 
derived using the adopted laboratory wavelengths. Ob- 
servation al errors for the line bisectors were computed 
following [Gra3 (|1983f l. 

The observed line bisectors are shown in Fig. [TJ where 
we emphasize the location of the line cores. The cores of 
the strongest lines are blucshiftcd by about —100 ms" 1 
while those of weaker lines show more blueshift, up to 
—900 ms" 1 (note that the absolute velocity scale was 
determined using 3D model predictions, as explained in 
Sect. H|) . The overall shape of the bisectors resembles the 
letter C, although it varies with line strength. Strong 
(weak) lines resemble more the lower (upper) part of 
the letter C. This is not very different from the case of 
cool metal-rich stars, implying that low metallicity does 
not affect the ba sic behavior of bisector shapes (see also 
I Gray e t al. 2008|). The strongest lines in the spectrum of 
HD 122563 (residual flux of about 0.2) are more symmet- 
ric and their bisectors have a span of about 100 ms" 1 . 
The weaker lines (residual flux of about 0.8), on the other 
hand, have a bisector span of about 800 ms" 1 . In a hy- 
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drostatic atmosphere, the lines would be symmetric and 
not shifted. This hypothetical case is illustrated with the 
dashed line in Fig. Q] 

3. HYDRODYNAMIC SIMULATION 



Collet et al. (in preparat ion; see alsolCollet et al. 



2009) 



used the stagger code bv lNordlund fc Galsgaardl (|1995[ ) 
to compute a 3D radiative-hydrodynamic model with 
stellar parameters T c s = 4600 K, log g = 1.6 [cgs], and 
scaled solar chemical composition w ith [Fe/H] = —3. The 
stellar parameters were taken from lBarbuv et all ()2003[ ) 
with a slight revision of the surface gravity based on the 
parallax determin ation of HD122563 f rom the Hipparcos 
data reduction by Ivan Leeuwenl (|2007fl . 

The 3D model was computed using a discrete 
480x480x240 grid representing a cube of physical dimen- 
sions 3700x3700x1100 Mm 3 . The equations of mass, en- 
ergy, and momentum conservation were solved together 
with the radiative transfer equations adopting realistic 
equation of state and opacities. The simulation sequence 
used for the present work consists of 51 snapshots that 
cover about 150 hours of stellar time and it typically 
contains about 10 granules at the surface. The evolution 
timescale o f the granul e s is a bout 90 hours. For more 
details, see lCollet et~aTI (|2009D . 

We used the 3D model to synthesize all 58 Fel lines 
shown in Fig Q] For the line synthesis we adopted a 
coarser grid of 50x50x125 points to which we interpo- 
lated the original simulation. A number of tests were 
performed with a higher-resolution grid of 100x100x125 
points to ensure that this had no significant impact on 
the shapes of the synthetic line profiles. 

4. LINE PROFILES: SYNTHETIC VS. OBSERVED 

The synthetic 3D profiles were convolved with a Gaus- 
sian instrumental profile of FWHM = 1.5 kms -1 , which 
corresponds to the spectral resolution of our data. From 
a x 2 line profile fitting of about 30 Fe I lines with the high- 
est local S/N, with Vsini and iron abundance as free 
parameters, we obtained Vsini = 3.2 ± 0.6 kms -1 . For 
the y 2 fitting w e convolved the disk-integrated profiles 
(e.g.. lGravll992l p. 370). To calculate the disk-integrated 
profiles for comparison with the observational data we 
applied a Vs'mi — 3.2 kms -1 to the emergent intensi- 
ties. A detailed abundance analysis of our HD 122563 
spectrum using the 3D model will be provided in a later 
publication^ 

In Fig. [2] we show a comparison of observations and 
model for 3 representative Fel lines. The observed lines 
have been shifted by about +30, +30, and -100 ms -1 
(from strongest to weakest line) to make their core ve- 
locities agree with the theoretical values. This allows 
a better comparison of line bisector shapes because it 
removes the uncertainties in the measured core wave- 
lengths and rest laboratory wavelengths. We discuss the 
core wavelength shifts below and show that the artificial 
shifts applied here are within the observational errors Q 
The average difference in bisector velocities (model minus 

2 According to the preliminary analysis bv lCollet et al. (2009), 
the 3D LTE iron abundances determined from Fel lines show a 
trend with excitation potential, going from about loge(Fe) = 4.0 
at x = OeV to loge (Fe) = 4.6 at \ = 4eV, with some line-to-line 
scatter comparable to the one from classical ID analyses. 

3 We avoid the comparison of line bisectors averaged in groups 
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Figure 2. Bisectors of three representative Fel lines in the spec- 
trum of HD 122563 (solid lines). Error bars (~ 30 — 150 ms from 
core to continuum) have been omitted for clarity. The dot-dashed 
lines correspond to 3D model predictions. The observed lines in 
this figure have been shifted to match the theoretical line cores and 
thus allow a better comparison of bisector shapes. The dashed line 
at zero velocity represents ID model predictions. 

observations for all 58 Fel lines) is 10 ± 71 ms -1 (here- 
after, error bars correspond to the 1 a scatter). This 
number, however, depends on the relative flux. Near 
the continuum, the bisector point-by-point difference is 
about —30 ± 145 ms -1 while it is 9 ± 33 ms" 1 at a 
residual flux of about 0.2. At all flux values the differ- 
ence is consistent with zero within the 1 a scatter. The 
formal error in the observed bisector velocities shows a 
residual flux dependency remarkably similar to that de- 
scribed above; i.e., it ranges from about 30 to 150 ms -1 
from core to continuum. Thus, the scatter seen in the 
differences between theoretical and observed bisectors is 
fully compatible with the 1 a observational errors. 

Measurements of core wavelength shifts (Av c ) are in 
principle more precise than line bisectors, since only a 
small portion of the line profile is examined, thus lower- 
ing the probability of being affected by blends. In Fig. [3] 
we show the relation between Av c and line strength . The 
so-ca lled "third signature" of stellar granulation (jGravi 
120091 ). i.e., the increased blueshift with decreasing line 
strength, is clearly detected in the data and very closely 
reproduced by the 3D model. The third signature is 
due to the fact that weaker lines form in deeper layers, 
where the granulation velocities and intensity contrast 
are larger. The line-to- line scatter seen in Fig. [3] increases 
towards weaker lines, an effect that is present also in the 
3D model predictions. It is tempting to conclude that 
this increased scatter is real, but given the timescale of 
the evolution of granules in comparison with the dura- 
tion of the simulation and the possibility that line cores 
of weaker lines are more difficult to measure (but see be- 
low), we must refrain from making a strong statement 
about this finding. Note also that the variation of the 
continuous opacity over the wavelength range covered 
by our data is substantial, which introduces wavelength 
dependent sh ifts that might contr ibute to the observed 
scatter (e.g.. iDravins et aLl 119811) . The difference be- 
tween observations and model predictions for Av c has 
alcr scatter of 104 ms" 1 . 

of line strength because this leads to a certain loss of information. 
The detailed shapes of lines of similar strength can be different due 
to differences in the properties of the atomic transition involved. 
Averaging is justified in cases where the S/N is low, the spectral 
lines are s everely blended, or t he granulation signatures are too 
weak (e.g., Ramirez et al. 2008). 
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Figure 3. Core wavelength shifts of 58 Fel lines in the spec- 
trum of HD 122563 (filled circles) and those predicted by the 3D 
model (open squares) as a function of residual flux (top panel) and 
equivalent width (bottom panel) . The residuals (observation minus 
model) are shown with crosses. The dashed line is at zero velocity. 

To estimate the error in the measured Av c values we 
artificially introduced Gaussian noise to the theoretical 
profiles to simulate the local S/N of the data. We took 
into account the variation of the S/N across the line 
profile since it varies according to the amount of flux de- 
tected at each wavelength and it is therefore lower at the 
line core compared to the local continuum. This experi- 
ment was performed after degrading the theoretical pro- 
files to the spectral resolution of our data and convolving 
it with a rotational profile with the V sin i we derived for 
HD 122563. Finally, for each line tested we added in 
quadrature the error introduced by the uncertainty in 
the measured laboratory wavelengths (20 ms -1 ) and the 
estimated error of the wavelength calibration of the data 
(30ms -1 ). The Av c values obtained from 1,000 of these 
tests, averaged over the 58 Fel lines, have a 1 a scatter 
of ~ 108 ms -1 , nearly independent of line strength but 
mildly dependent on wavelength, reflecting the variation 
of the S/N across our observed spectrum. This estimated 
observational error is very similar to the scatter seen in 
the comparison of model predictions and observations. 

5. SUMMARY AND DISCUSSION 

We have observed the very metal-poor red giant star 
HD 122563, a primer for studies of the Galactic stellar 
halo, at very high resolution and high signal-to-noise ra- 
tio. The high quality of our data allowed us to measure 
the asymmetries and core wavelength shifts of a num- 



ber of Fel lines, which clearly revealed the signatures 
of granulation. A state-of-the-art 3D LTE radiative- 
hydrodynamic simulation was used to synthesize the 
same set of lines. Using line profile fitting we derived 
a projected rotational velocity V svai = 3.2 ± 0.6 kms -1 
while the analysis of core wavelength shifts allowed us 
to infer an absolute radial velocity V r — —25.39 ± 
0.09 kms -1 . The line bisectors and core wavelength 
shifts computed from the 3D model predictions agree 
with the observations of HD 122563 remarkably well. The 
average difference between observations and model is 
consistent with zero within the estimated 1 a observa- 
tional errors. This suggests that the 3D model contains 
realistic information regarding the atmospheric structure 
and velocity fields of HD 122563. 

Abundance analyses of HD 122563 using our 3D model 
show that standard ID estimates of elemental abun- 
dances caii_be_off try up to one order of magni- 
tude (jCollet et all 120071 120091) , Qualitatively simi- 
lar results have been obtained with a simulation com- 
puted by an independen t group (jlvanauskas et alj 120101 : 
IDobrovolskas et aT1l2010f ) , although their abundance cor- 
rections tend to be smaller. The large 3D corrections will 
have an impact on our interpretations of stellar abun- 
dance data for halo stars and may lead us to revise their 
connection with early enrichment of the Galaxy, super- 
novae yields, and internal mixing processes. The observa- 
tions presented in this Letter and their comparison to the 
3D model predictions suggest that such large 3D effects 
are reliable and must therefore be seriously considered in 
all studies of metal-poor stellar abundances. 

Admittedly, the largest 3D corrections apply to abun- 
dances derived from molecular features whereas we have 
analyzed only Fel lines. The Fel lines provide one 
sampling of the atmospheric structure with a particu- 
lar weighting over the different structures and velocity 
fields. Other species, particularly the molecules, will 
provide different weighting and should in principle be 
investigated before claiming that the 3D models are ad- 
equately justified for all abundance analyses. Unfortu- 
nately, in our spectral range the molecular lines are too 
weak for this type of work. We note, however, that our 
Fe I lines form over a wide range of depths, spanning a 
large portion of the photosphere, including the regions of 
formation of molecular lines typically used in abundance 
work (logTR OSS ~ -2). 

Both our 3D model and line synthesis assumed LTE, 
yet the impact of departures from the LTE approxima- 
tion is likely to be very important at low metallicities 
(e.g.. IAsplundl[2005l) . To first order, however, non-LTE 
affects the abundances but not the shapes of line pro- 
files. A fully consistent 3D non-LTE model atmosphere 
and line formation scheme is currently beyond our reach. 
Nevertheless, our work represents one step forward in the 
direction of major improvements to elemental abundance 
determination techniques. 

This work was supported in part by the Robert A. 
Welch Foundation of Houston, Texas (grant number F- 
634). We thank the anonymous referee for helping us 
improve the original manuscript. 
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